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ABSTRACT 



Af-body simulations predict that dark matter haloes are described by specific density profiles on both galactic- and cluster-sized 
scales. Weak gravitational lensing through the measurements of their first and second order properties, shear and flexion, is a powerful 
observational tool for investigating the true shape of these profiles. One of the three-parameter density profiles recently favoured in 
the description of dark matter haloes is the Einasto profile. We present exact expressions for the shear and the first and second flexions 
of Einasto dark matter haloes derived using a Mellin-transform formalism in terms of the Fox H and Meijer G functions, that are valid 
for general values of the Einasto index. The resulting expressions can be written as series expansions that permit us to investigate the 
asymptotic behaviour of these quantities. Moreover, we compare the shear and flexion of the Einasto profile with those of different 
mass profiles including the singular isothermal sphere, the Navarro-Frenk- White profile, and the Sersic profile. We investigate the 
concentration and index dependences of the Einasto profile, finding that the shear and second flexion could be used to determine the 
halo concentration, whilst for the Einasto index the shear and first and second flexions may be employed. We also provide simplified 
expressions for the weak lensing properties and other lensing quantities in terms of the generalized hypergeometric function. 
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1. Introduction 



A more accurate description of the elements that constitute 
our universe, such as the dark matter haloes that are believed 
to exist around galaxies and clusters, is of crucial importance 
for our understanding of cosmological structural formation. 
Recent results from A^-body simulations of cold dark matter 
(CDM) (Navarro et al. 2004; Merritt et al. 2006; Gao et al. 2008; 
Hayashi & White 2008; Stadel et al. 2009; Navarro et al. 2010; 
Reed et al. 201 1; Vera-Ciro et al. 2012) indicate that nonsingu- 
lar three-parameter models such as the Einasto (1965) profile, 
fit a wide range of dark matter haloes better than singular two- 
parameter models, e.g. the Navarro, Frenk, and White (NEW) 
profile (Navarro, Erenk, & White 1996, 1997). 
The Einasto profile is given by 



p(r) = ps exp < 











-'1 


1 









(1) 



where r is the spatial radius, the shape parameter n is called the 
Einasto index, represents the radius of the sphere that con- 
tains half of the total mass, ps is the mass density at r = rj, 
and d„ is a function that ensures that is indeed the half- 
mass radius. An analytical expansion for the function d„ » 
3n- l/3-i-8/1215n-i-<9(n^) is provided by Retana-Montenegro et 
al. (Retana-Montenegro et al. 2012, hereafter RM12). One im- 
portant characteristic of this profile is that its power-law loga- 
rithmic slope, yir) = -dlnp/dlnr ~ r'^", depends on the Einasto 
index, which provides a profile that more accurately fits in the 
inner regions of simulated dark matter haloes than other profiles 
such as the NEW profile. In the study of real galaxies, several 
authors have used multi-component Einasto models, consisting 
generally of two or more Einasto components for each galaxy. 



where each component represents a homogeneous stellar pop- 
ulation with its own set of parameters. Eor example, some of 
the first galaxies to be modelled using multi-component Einasto 
models were M31 by Einasto (1969a) with values of 0.25 < 
n < 1, and other nearby galaxies such as Milky Way, M87, 
M32, Eornax, and Sculptor, and M31 by Einasto (1974) with 
0.5 < « < 4. Later, in a series of papers multi-component Einasto 
models were employed to model the luminous components of 
several galaxies such as the Milky Way (Einasto & Haud 1989), 
M87 (Tenjesetal. 1991), M31 (Tenjes et al. 1994), and M81 
(Tenjes et al. 1998); in these papers, the Einasto index is char- 
acterised by values of 0.36 < n < 7.1. The seven distant 
spiral galaxies GSS 074-2237, GSS 064-4412, GSS 094-2210, 
GSS 104-4024, GSS 064-4442, MDS uemO-043, and HDFS 
J223247.66-603335.9 were studied by Tamm & Tenjes (2003) 
and Tamm & Tenjes (2005), respectively. As in the earlier works 
mentioned, they modelled each visual component with a Einasto 
profile, the authors found values of 0.25 < n < 0.91 and noted 
that the Einasto indices for the disk component of the galax- 
ies at high redshift follow a trend of having smaller values 
than the ones at lower redshift. Tempel & Tenjes (2006) fitted 
a multi-component Einasto model to the Sombrero galaxy, with 
0.78 < « < 3 for the visual components. Tamm et al. (2007) and 
Tempel et al. (2007) presented a multi-component Einasto law 
study of M3 1 : using photometric data and metallicity measure- 
ments, they obtained the matter distribution of luminous com- 
ponents with 0.70 < n < 4.20, then tried to fit several mod- 
els for the dark matter halo using kinematical data from the 
literature to construct a dynamical model and derive the dark 
matter density of the galaxy, they concluded that Einasto and 
NEW profiles give the best fits. Dhar & Williams (2011) fitted 
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the surface brightness density of a sample of elHptical galax- 
ies using a multi-component Einasto profile, finding values of 
1 < n < 3 for the central components, and 5 < n < 8 for the 
outer components. Chemin et al. (2011), who studied the rota- 
tion curves of low mass spiral galaxies, modelled the dark mat- 
ter halo with a Einasto profile, and obtained smaller values of n 
than predicted by computational simulations. On the other hand, 
according to A^-body numerical calculations the Einasto index 
dependens on both the halo mass and redshift (Hayashi & White 
2008; Gao et al. 2008). Typical values of the Einasto index are 
in the range 5 < n < 8 according to the results of A^-body sim- 
ulations (Navarro et al. 2004; Gao et al. 2008; Hayashi & White 
2008; Navarro et al. 2010). Vera-Ciro et al. (2012) analysed dark 
matter haloes of Milky Way-like systems and concluded that the 
Einasto model with values of 2 < n < 5 is preferred over the 
NEW profile. 

An alternative form of the density often used in dark matter 
halo studies is 



p (r) = p_2 exp \ -In 



l/n 

— I -1 

r-2 



(2) 



where r_2 is the radius at which the logarithmic slope of the den- 
sity distribution has a value of -2 and p_2 - pir-i)- A useful 
quantity to define is the concentration ce = ?"20()/'"-2, where r2oo 
is the virial radius of a halo of mass M200, whose density is 200 
times the critical density of the Universe at the halo redshift. 
One of the advantages of the Einasto profile over other profiles 
is that it has excellent agreement with the conditions outlined by 
Einasto (1969b) for constructing real galactic models, specifi- 
cally, some moments must be finite. In particular, for this profile 
some moments, such as the total mass, central gravitational po- 
tential, and effective radius, are finite. In contrast, other profiles 
have logarithmic moments that must be truncated at some radius 
to ensure that the profile remains finite. 

Gravitational lensing provides a direct way to study the 
mass distribution of large structures in the universe, such as 
galaxies and clusters, without making any assumptions about 
their dynamical state or composition. Lensing studies taking 
advantage of high-quality imaging have proven to be suc- 
cessful in mapping the distribution of dark matter in clusters 
and galaxies (Kneib et al. 2003; Clowe et al. 2006; Bradac et al. 
2008; Abate et al. 2009; Okabe et al. 2010; Oguri et al. 2010; 
Okabeetal. 2011; Huang et al. 2011; Ragozzine et al. 2012; 
Jee et al. 2012; Oguri et al. 2012). There are two lensing 
regimes: the strong regime, where multiple images or strong 
distortions of a galaxy can be produced by an interven- 
ing distribution of matter, and the weak regime, where the 
lensed galaxy image is only slightly distorted, causing the 
intrinsic elliptical galaxy to appear as a distorted elliptical 
image. Weak lensing is a valuable and accurate tool for 
determining the shapes of dark-matter-halo density profiles, 
such as ellipticity (Hoekstra et al. 2004; Mandelbaum et al. 
2006; Parker etal. 2007; Evans & Bridle 2009; Deb et al. 
2010; Oguri etal. 2010) and triaxiality (Oguri et al. 2005; 
Gavazzi 2005; Corless&King 2007; Corless et al. 2009; 
Feroz & Hobson 2012). Until now, most weak lensing studies 
have considered only linear-order effects such as the weak shear, 
the quantity responsible for the induced ellipticity in the galaxy 
(see e.g. Kaiser et al. (1995); Bartelmann & Schneider (2001); 
Schneider et al. (2006); Hoekstra & Jain (2008) for reviews). 

In the past few years, the study of high-order lensing prop- 
erties has grown in importance (Goldberg & Natarajan 2002, 
Goldberg & Bacon 2005, Bacon et al. (2006, hereafter B06)). 
These properties written as high-order derivatives of the deflec- 



tion potential can be recognized as convergence and shear gra- 
dients. The convergence gradient, called the first flexion 'F, in- 
duces a centroid shift in the lensed image with respect to the 
source or "skewness". The shear gradient, called second flex- 
ion 0, generates an arc-like shape in the lensed image or "arc- 
ness". Weak flexion provides useful information about dark mat- 
ter haloes on galactic- and cluster-sized scales, particularly when 
probing substructure on smaller-scales where flexion is more 
sensitive to shear-only studies (Leonard et al. 2009; Bacon et al. 
2010; Er etal. 2010). 

Several methods have been developed to measure the flex- 
ion of a lensed image, e.g. shapelets (B06; Massey et al. 
2007; Fluke et al. 2012) and surface brightness moments 
(Irwin & Shmakova 2006; Irwin & Shmakova 2005; Irwin et al. 
2007; Goldberg & Leonard 2007; Okuraetal. 2007, 2008; 
Okura & Futamase 2009; Schneider & Er 2008). Cain et al. 
(2011) introduced a new method, called the analytic im- 
age model (AIM), to study flexion in astronomical im- 
ages. Observational measurements of flexion include the de- 
tection of mass substructure in the Abell 1689 cluster by 
Leonard etal. (2007) and Cain et al. (2011) using observa- 
tions of the Hubble Space Telescope (HST), as well as 
Okura et al. (2008) employing Subaru images; galaxy-galaxy 
flexion detection in the ground-based survey Deep Lens Survey 
(Goldberg & Bacon 2005) and the space-based HST COSMOS 
survey (Velander et al. 201 1). 

In addition, flexion has been proposed as a powerful cosmo- 
logical tool: Munshi et al. (2011) suggested the use of conver- 
gence shear and flexion maps to decrease errors in the measur- 
ing standard candles distances. Camera & Diaferio (2011) stud- 
ied how the flexion signal-to-noise ratio could be used to dis- 
cern between cosmological models, Munshi et al. (2011) and 
Schafer et al. (2012) proposed the use of cosmic flexion to probe 
large-scale structure. Hawken & Bridle (2009) studied the halo 
ellipticity on galactic scales, and found that the inclusion of flex- 
ion yields tighter constraints on ellipticity than shear-only stud- 
ies. Er & Schneider (201 1) and Er et al. (201 1) proposed a new 
way to determine the halo ellipticity using the ratio of tangential- 
to-radial flexion and studied its behaviour as a radius function. 
Er et al. (2012) concluded that flexion is more sensitive to ellip- 
ticity than shear by performing a likehood analysis of mock flex- 
ion and shear data. Additionally, Viola et al. (2012) considered 
the case in which cross-terms between both shear and flexion 
and between intrinsic galaxy ellipticities and flexion are not ig- 
nored, concluding that these terms can cause a considerable bias 
in the flexion estimations. 

In view of the increased use of the Einasto profile 
in cosmological studies (e.g. see Catena & Ullio (2010); 
Reed etal. (2011); Chemin et al. (2011); Dhar & Williams 
(201 1); Catena & Uflio (201 1); Narikawa & Yamamoto (2012)), 
it is natural to extend its applications to weak shear and flex- 
ion lensing studies. Previously, several authors had performed 
weak lensing studies using the Einasto profile. For example, 
Hayashi & White (2008) measured the cross-correlations be- 
tween halo centres and mass, and between galaxies and mass, 
in the Millennium Run (Springel et al. 2005), and found that 
the Einasto profile provides a close fit in the inner regions of 
their two-part model of the halo-mass cross-correlation func- 
tion. Mandelbaum et al. (2008) analysed, using a weak statis- 
tical approach, a sample of galactic- and cluster-sized dark mat- 
ter haloes from the Sloan Digital Sky Survey, and obtained very 
similar concentration-mass relations for the NFW and Einasto 
profiles. Mamon et al. (2010) used analytical approximations of 
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the shear of the Einasto profile to compare it with the NFW 
shear. 

Parametric models such as the singular isothermal sphere, 
the NFW, and Sersic profiles have been used to model 
the dark matter distribution in weak lensing analyses (e.g. 
Hawken & Bridle 2009; Umetsu et al. 2011; van Uitert et al. 
201 1; Er et al. 201 1; Viola et al. 2012; Soucail 2012), the prop- 
erties of these models having been studied by several authors 
(Wright & Brainerd 2000; B06; Lasky & Fluke 2009). In the 
case of the Einasto profile, RM12 studied the analytical prop- 
erties of the Einasto profile by applying a Mellin-transform for- 
malism. In terms of Fox H and Meijer G functions, they derived 
analytical expressions of lensing properties for all values of the 
Einasto index, concentrating on the surface mass density, cumu- 
lative mass, deflection angle, and deflection potential. However, 
by means of the Mellin-transform formalism it is possible to 
extensively study the weak-lensing analytical properties of the 
Einasto profile. This study provides analytical expressions that 
can used to model realistic Einasto dark matter haloes in weak 
lensing modelling studies. 

In this work, we apply Mellin-transform formalism to ob- 
tain and study in detail the analytical expressions for the weak 
lensing properties of the Einasto profile: the shear, and first, and 
second flexions. This paper is organized as follows. We sum- 
marize the weak lensing formalism in Section 2.1, and present 
the Mellin-transform technique in Section 2.2. In Section 3 we 
derive closed expressions for the shear and first and second flex- 
ions in terms of the Fox H and Meijer G functions. We then 
use the series expansions of these expressions to investigate their 
asymptotic behaviour In Section 4, we compare our results with 
those for the SIS, NFW, and Sersic profiles. In Section 5, we 
summarise and discuss our main results. Finally, in the appen- 
dices A and B we provide series expansions of the lensing prop- 
erties, and explicit expressions in terms of the generalized hy- 
pergeometric function, respectively. Throughout the paper, we 
adopt a cosmological model with the matter density Q.m - 0.26, 
the cosmological constant Qa - 0.74, and the Hubble constant 
Hq = 72kms-'Mpc"'. 



2. Theory 

We provide a brief description of the two main theoretical as- 
pects employed throughout this paper. 



2.1. Weak lensing formalism 

The weak lensing formalism using complex notation was intro- 
duced by B06. In the thin lens approximation, the lens equation 
is given by (Schneider et al. 1992) 



jg = 6» - Viff(6), 



(3) 



where j8 and denote the positions on the source plane, and on 
the image plane, respectively, and iff{ff) is the deflection potential 
defined by a two-dimensional Poisson, V^iff(ff) - 2k{9), with the 
convergence k(0). Moreover, the convergence can be written as 



where 1,(6) is the surface mass density. 



47rGDLZ)Ls 



(4) 



(5) 



is the critical surface mass density, and Dl, Ds, and Dls are the 
angular distances from observer to lens, from observer to source, 
and from lens to source, respectively. In addition, it is convenient 
to define the complex gradient operator (B06) 

5 = — + ! — = de'f, (6) 
56*1 802 

where is the rotation angle, relative to the basis. The d operator 
is simply a spin-i' raising operator and its complex conjugate d* 
a spin- 5 lowering operator. 

When we study gravitational lensing on scales where the deflec- 
tion potential changes are larger than the scale of the lensed im- 
age, we can expand up to second order the lens equation in eq. 
(3) around the neighbourhood of the lensed image^ 



/3i = ^HijOj + ^Dijkejek, 

where is the Jacobian matrix defined by 



l-K-yi -72 



-72 1 - /c + 71 / ' 

with the convergence. 



(7) 



(8) 



(9) 



which is a spin-0 field, and with Ji - \ ((Ax.i + ^\y),y2 - the 
components of the complex shear 



r = n + ' n = Irl e^"* = ^<95iA, 



which is a spin-2 field. The matrix Dijk - dAijIdOk describes the 
behaviour of the convergence and shear across the lensed image 
by introducing two new lensing properties 



i^ijk - Tijk + Gijk, 



(11) 



namely the first flexion field or spin-1 first flexion Tijk and the 
second flexion field or spin-3 second flexion Qijk. Both field 
components can be expressed as third-order derivatives of the 
deflexion potential (Hawken & Bridle 2009) 



Q2^\ [^^xxy - lA.vvv) . 



(12) 

(13) 

(14) 

(15) 



and taking advantage of the complex formalism, we can com- 
pactly write the first and second flexions as 

r = ri + ir2^^ddd*ifr, (i6) 



+ 102 = ]j^ddd<p. 



(17) 



From eqs. (16) and (17), one can clearly see the rotation symme- 
try for both flexions. Applying the complex conjugate operator 



' We define a spin-^ lensing quantity by requiring that it is invariant 
under rotations = Inj s, where s is any natural number except zero. 
- We do not consider crossed terms in the lens equation expansion. 
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to eqs. (9), and (10) and comparing with eqs. (16), and (17), we 
find a compact and elegant definition of the second-order prop- 
erties as gradients of the first-order lensing properties 

r ^\r\e''l' ^dK^d*y, (18) 



(19) 



2.2. Mellin-transform technique 



3.1. Convergence and shear 

In the weak lensing regime up to first order, the lensed galaxy 
image has two distortions: the convergence k causes an isotropic 
stretching in the lensed image, which magnifies the image by 
increasing its size, and the shear y also causes an anisotropic 
stretching in the lensed image, that is responsible for the induced 
ellipticity in the lensed galaxy. 

Foremost, to calculate the convergence, we must project the den- 
sity profile on the lens plane using an Abel transform 



The Mellin transform technique (Marichev 1983; Adamchick 
1996; Fikioris 2007) consists in that one-dimensional definite 
integrals 



2(^)=2 



r 



p(r) rdr 



(25) 



fiz) 



Jo 



g{t,z)dt. 



(20) 



where ^ is the radius from the lens centre, and r is the spatial 
radius. We follow here the notation of RM12 for the Einasto 
profile 



can be expressed as the Mellin convolution of the functions /i 
and/2 



p(r) = po exp 



fiz) 



where we define the central density po - Ps e''" - p-i e 
scale length h - r^ld"„ - r-2l (2n)". Additionally, we define 
another quantity, the central convergence 



(26) 



2" and 



The Mellin convolution theorem, which states that the Mellin 
transform of a Mellin convolution of two functions is the point- 
wise product of their Mellin transforms, can be applied to eq. 
(21) inverting the Mellin transform of the Mellin convolution, 
f(z) can be expressed as the inverse Mellin transform of the 
pointwise product of the f\ and /2 Mellin transforms. The Mellin 
transform is defined by 



E(0) 2pQhnr(n) 



Combining eqs. (4), (25), (26) and (27), we get 



k(x) 



r'O. 

nr(n) X 



sds 



(27) 



(28) 



r 

u) - 4>(u) - I 
Jo 



f(z)z"-'dz, 



(22) 



and the inverse Mellin transform by 



where x - 6D]^/h - ^/h and s - rjh are the dimensionless radii. 
The integral in eq. (28) cannot be expressed in terms of elemen- 
tary or special functions for general values of n. However, using 
the Mellin transform technique explained in Section (2.2), we 
can write this integral as a Mellin-Barnes integral 



m,\z)^m^^.Uu)z-^Au, (23) ^^^^ J rTi2ny)T{-\^y) 

k(x)^ —J 



where the integration path is a vertical line in the complex plane. 
The integral in eq. (20) may then be written as 

fiz) = ^ J^9JV,(«)9%(")^"'dM- (24) 

With the requirement that f\ and fo are of hypergeometric type, 
their Mellin transforms can be written as products of the form 
T{a+ Au) or [F (a + Am)]"', with F (v) the gamma function and 
A real. The resulting integral in eq. (24) is of the Mellin-Barnes 
type and it then can be evaluated as either a Fox H function for 
A 1 or as a Meijer G function for A = 1 . 



3. Weak lensing 

We derive closed expressions for the weak-lensing first- and 
second-order properties of the Einasto profile: the shear y, and 
the first T and second Q flexions in terms of Fox H, Meijer G 
functions, and the generalized hypergeometric function. Using 
these expressions, we calculate the expansion series and investi- 
gate its asymptotic behaviour. The results of this section provide 
a useful and straightforward way to study weak lensing, where 
the matter distribution is believed to be described by an Einasto 
profile. 



F (n) 2ni 



r(y) 



[^Ydy. (29) 



The Fox H function (Fox 1961) is denoted as a Mellin-Barnes 
integral. 



Tjm,n 



(a, A) 




ib,B) 


z 



-f 



n"ir(^ + fi,..)n;=ir(i-fl,-A,..) 

n;i,„,, r(l - bj - Bjs) U%„,, riaj + Ajs) 



(30) 



Comparing the integral in eq. (29) with the above definition, we 
obtain a close expression for the convergence in terms of the Fox 
if function (RM12) 



yfn 2,0 
K (x) - _ , . X H 



Tin) 



1,2 



(0, 1) 
(0, 2«),(-i, 1) 



(31) 



The shear for an circularly symmetric lens is (Miralda-Escude 
1991) 

y{x)-K (x) - k{x) , (32) 



where 



2 r 

k(x) = — x'K(x')dx', (33) 
Jo 
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is the average convergence within the dimensionless radius x. 
Inserting eq. (31) into eq. (33) and substituting this result along 
with eq. (31) into eq. (32), we can re-express the resulting in- 
tegral as a Fox H function and obtain the shear for the Einasto 
profile 



with 



2,1 



(-5, IX (0, 1) 
(0, 2n),(i, 1),(-|, 



1) 



(34) 



Eq. (34) provides an expression for the shear in terms of one Fox 
H function, instead of two Fox H functions as found by RM12. 
Writing eq. (34) in terms of one instead of two Fox H functions 
makes it easier to manipulate for analytical and numerical pur- 
poses. 

The Fox H function is a very general function and reduces to 
most of the elementary and special functions. Despite not being 
a common special function, it has great potential as an analytical 
tool in theoretical astrophysics, in particular to study the ana- 
lytical properties of density models such as the Sersic profile 
(Baes & Gentile 201 1; Baes & van Hese 201 1) and Einasto pro- 
file (RM12). This function will be included in future versions of 
the software Mathematica. Additionally, several authors such 
as Yilmaz & Alouini (2009) and Shafique Ansari et al. (2012) 
have made available accurate and fast numerical routines to com- 
pute the Fox H function. Details about the many properties of 
this function can be found in Mathai (1978), Srivastava et al. 
(1982), Kilbas & Saigo (2004), and Mathai et al. (2009). 
We remark that the shear y(x) is not a directly measur- 
able property owing to the mass-sheet degeneracy, but that 
the measurable property is the reduced shear (Falco et al. 
1985; Gorenstein et al. 1988; Schneider & Seitz 1995; 
Narayan & Bartelmann 1996) 

gix)^-, TT- (35) 

1 - k{x) 

3.2. First and second flexions 

Considering the weak lensing regime up to second order, two 
new lensing properties can be recognized: the first flexion f, 
which describes the behaviour of the convergence gradient 
across the lensed image and the second flexion ff, which de- 
scribes the behaviour of the shear gradient across the lensed im- 
age. These flexions produce centre-shift and arc-like distortions 
that, with the addition of the shear, cause the lensed image of the 
elliptical galactic source appear to have a "banana-like" shape. 
The first flexion can be found by simply calculating the conver- 
gence gradient, using eqs. (31) and (18) 



G' (x) = H: 



(-1,2),(-1, 1),(0, 1),(1,2) 
2n),(-i 1),(2, 2),(0, 2),(-|, 



1) 



(38) 

We have obtained analytical expressions for both flexions ex- 
pressed as Fox H functions, where the flexions are circularly 
symmetric as expected for the Einasto profile. The second flex- 
ion with nine gamma functions in the integrand is a more compli- 
cated function than the first flexion with only five gamma func- 
tions. 

As indicated before, the shear is affected by the so-called mass- 
sheet degeneracy, and the same measurement difficulty arises 
and with the first and second flexion. Schneider & Er (2008) 
demonstrated that the observable properties are the reduced flex- 
ions 



and 



F{x) 



G(x) = 



r(x) + g(x) r*(x) 

I-k(x) 

g(x) + g(x)g(x) 

I-k(x) 



(39) 



(40) 



where F (x) is a spin-0 field, G (x) is a spin-3 field, and 'F* (x) 
is the complex conjugate of the first flexion. 
Given that the scope of this work is to present analytical expres- 
sions of different weak lensing properties and the quantitative 
comparison of these properties for several profiles, we focus on 
the lensing properties given in eqs. (34), (36), and (37), instead 
of the reduced ones in eqs. (35), (39), and (40). 

3.3. Integer and half-integer values of n 

We simplify the expressions for the shear in eq. (34) and first, in 
eq. (36), and second, in eq. (37), flexions in terms of the Fox H 
function for rational values of «, to the case when n is an integer 
or half-integer number, where the resulting expressions can be 
written in terms of the Meijer G function. 
The Meijer G function is defined by the Mellin-Barnes integral 
(Meijer 1936) 



a 




b 


z 



- f- 

Jx n 



u"unbj + s) n;;=ir(i-fl,-i) 



9 

j=m+ 



■s)U%„,J(aj + s) 



z'-'ds. (41) 



r(x)^ro {H: 



'2,3 



with 



(-1,2),(0, 1) 
(0, 2«),(-i, 1),(0, 2) 



(36) 



hr(n) 



the flexion amplitude. 

Combining eqs. (19) and (34), plus some algebra, the second 
flexion may be obtained 



There is extensive literature about the Meijer G function and its 
many useful properties (Bateman & Erdelyi 1953; Luke 1976; 
Andrews 1985; Prudnikov et al. 1990). 

For the Fox H function, for there is not yet a numerical imple- 
mentation contrarily, there are various software packages with 
Meijer G numerical routines, such as the commercial Maple, 
Mathematica, and the free open-source Sage and mpmath li- 
brary. 

By substituting the Gauss multiplication formula 
(Abramowitz & Stegun 1970) 



g(x)^-^G'(x)e'"^, 



2n-\ 



(37) r(2ny) = (2n)-5+2«.v {2n)^-n Yiy) p] r + , (42) 
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into eq. (29), and making use of eq. (41) for the comparison, we 
have (RM12) 



k(x) 



k(x) 



-,2/7,0 



12: 

2(27r)"-i Vnr(n) 



where ft is a vector of size 2n given by 
, f 1 2 2n- 



, 2n 2n 



2n 





x2 1 


b 


(2«)2«J 


1 















1 

_ 2' 


1 

2 


4 



(43a) 



(43b) 



\r(x)\^^Gli 



I 1^3 
2' 2' 2 
_1 

1 _? 1 

2' 2' 2 



I^(X)| 



'2' 2 



ill 

"2' 2' 2' 



(48) 
(49) 
(50) 
(51) 



which is an expression for the convergence of the Einasto profile 
in terms of the Meijer G function. 
Now, substituting the convergence into eq. (33) and perform- We may write the above Meijer G functions in terms of Bessel 
ing the integration of Meijer G function (eq. 07.34.21.0003.01 functions 
at the Wolfram Functions Site^), and inserting the integral prod- 
uct along with the convergence in eq. (32), we may write the 
shear as 



K (x) = K^X Ki (x) 



y(x) 



■XX 



2(27r)"-i ^/H^(n) 

r^2n,\ 
*-^l,2«+l 



b'. 



(In) 



2n 



(44a) 



T(x). 



(2;r)"-i ^lnn 



■,2n,l 
^I,2n+1 



(2n) 



2n 



2 3 

l-^K2{x)-^K, (x) 



\r(x)\^ ^xKoix), 

h - K2(x)\ - X Ko(x) - 4 Ki (x) 



(52) 
(53) 

(54) 
, (55) 



where b' is a vector of size 2n given by 

^, _ M 2 2n - 1 1 
1 2n ' 2« ' ' 2n ' 2 J 

Calculating the gradient for convergence following eq. (18) and fall-off 
using the diff'erentiation properties of Meijer G function (eq. 
07.34.20.0005.01 on the Wolfram Functions Site), we find 



with Ky(z) the modified Bessel function of the second kind of 
(44b) order v. 
For n ■■ 



J, the Einasto mass distribution presents a Gaussian 



p (r) = po exp 



(56) 



Therefore, we set n = 5 in eqs. (43), (44), (45), and (46) to 
(45) obtain 



Applying again the differentiation properties of Meijer G func- 
tion to derive the gradient for shear according to eq. (19), the 
result may be written as 



Six) 



To 



(27r)«-i V"^ 



■,2n+l,2 
'"3,2n+3 



_i __ 

2' 2' 2 
^ 3 i _3 



1_ _[_ 

23 

2' 2' 2 



(2«) 



2n 



(46) 



K(X) = /Cc X Gq J 



\j(x)\ ^ K^xG 



1,2 



\rix)\^2roG\ 



1,1 



-2ToG 



2' 2 
_l 

2 

_i i 

2' 2 

_i _i i 

2' 2' 2 
i 3 1. _3 
'2' 2' 2' 2 



(57) 
(58) 
(59) 
(60) 



Additional simplifications in terms of generalized hypergeomet 

ric function for half-integer values of n can be found in Appendix elementary functions' 



We can equivalently express these Meijer G functions in terms 



3.4. Simple cases: n - 1 and « = 5 

For n - \, the density profile decreases exponentially from the 
system centre 



K {x) = /fc e , 
lr(x)| = Je-- (e--l-x2), 

-2 5^0^6""', 



p (r) = po exp 



(-(e 



(47) 



\S(x)\ 



x^ 

\T{X)\: 

4 To 
— T- e 



e--l-^(xH2) 



(61) 
(62) 

(63) 

(64) 



For the exponential case, the resulting weak lensing expressions The results for the exponential profile described in eqs. (52), 
can be found by substituting n = 1 in eqs. (43), (44), (45), and (53), (54), and (55), and for the Gaussian profile in eqs. (61), 



(46) 



http://functions.wolfram.com/HypergeometricFunctions/MeijerG/ 



(62), (63), and (64) can be checked separately by substituting 
the density profiles in eqs. (47) and (56) into eq. (28) and per- 
forming the relevant derivations for the shear and flexions. 
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3.5. Asymptotic behaviour 

The behaviour of the weak lensing properties of the Einasto pro- 
file at small radii (x <K 1) can be studied by using the series 
expansions of Appendix B. We easily arrive at appropriate ex- 
pressions, which depend on the value of n for the shear 



r(i -») 2 
IrWI ~ -r — — -X 



IrWI 



4 r(l +n) 

r(3f) 



In On + 1) 



r(«)r(^) 



forn < 1, (65) 

forn = 1, (66) 
for/1 > 1. (67) 



The first flexion near zero behaves as 



\r(x)\ 



\r{x) 



\r(x)\ 



To 
n y/n 

To 

— r'' 
n yjT 

2«2 



r(l -n) X 



lng)+l+r 



r(^) 



forn<l, (68) 

for« = l, (69) 
forn>l. (70) 



For the second flexion, the asymptotic behaviour in the neigh- 
bourhood of the lensed image origin is described by 



\&(x)\ 



^° r(l-3n)x^ 



6n yfjT 



1 

for n < — or n 
3 



To 3 
y/jT 



for n 



1 

2' 
(71) 



1 

3' 
(72) 



1^(^)1 ~- 



35^0 



32« yfjT 



r(l -5n) x^ 



1 



2n2 \ 3n -H 1 



|^(x)|~^^[ln 
2n^ \ 3« -I- 1 



r(-^) 



x^l" forn>-, 
3 

(73) 

forn = 1, (74) 
for« > 1. (75) 



The power-law and logarithmic series converge very slowly at 
large radii {x » 1). Hence, we cannot use them to investigate 
the behaviour of the shear and flexions at large radii. However, 
following Kilbas & Saigo (1999), we derive asymptotic expan- 
sions at large radii for these properties. For the shear, we have 



,r(3«) _2 

2 -r: — — Kr- X 



2n 



Tin) Vnr(n) 
When jc — > oo, the first flexion behaves as 



-x"" 1--!- 
KcC X 2„ 



\r(x)\ 



„3/2 



and the second flexion is characterised by the behaviour 



\&(x) 



,3/2 



-L 

- e jc 2" . 



(76) 



(77) 



(78) 



4. Profile comparisons 

We compare the weak lensing properties for the Einasto pro- 
file obtained in the Section 3 with the properties of the 
singular isothermal sphere, the Navarro-Frenk- White profile, 
and Sersic model. The properties for these models can be 
found in Bartelmann 1996, Wright & Brainerd 2000, B06, and 
Lasky & Fluke 2009. We follow the same approach that was 
used by Wright & Brainerd (2000) and Lasky & Fluke (2009), 
which consists in fixing the halo mass M20() and permits the cal- 
culation of the virial radius. We use M2()o = 1 X 10i2/;-iMa 
and zi - 0.4, which are the approximate average galactic-sized 
halo mass and lens redshift, respectively, found by Parker et al. 
(2007) in their galaxy-galaxy weak lensing analysis of the 
Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). 
We choose to place the source at Zs - 0.92, which implies 
Dls/Ds - 0.5. We assume that the concentration-mass varies 
with the halo redshift in the range < z < 2 as (Duffy et al. 
2008) 



C200 (M20Q,Z) = A 



''^^(l+z)^ 



M, 



(79) 



where Mp, 



pivot 



2 X 10'2/z-iMo and I A, B, C] 



{6.40, -0.108, -0.62 ) for flie Einasto profile. Using the above 
relation, we find a concentration of ce - 5.80. In the case of the 
Einasto halo, we alse choose to use a value of n ^ 6 that cor- 
responds to the halo mass that is used according to Gao et al. 
(2008). For the NEW halo, the concentration can be deter- 
mined using eq. (79) with { A, B, C) = {5.71,-0.048,-0.47) 
(Duffy et al. 2008), we obtain cnfw = 5.31. The calculation of 
the parameters model is explicit, except in the case of the Sersic 
profile, for which we construct the model using the procedure 
outlined in Appendix B of Lasky & Fluke (2009). This proce- 
dure requires employing the Prugniel & Simien (1997) model, 
which is an analytical approximation of the Sersic deprojection; 
we use a Sersic index of m ^ 8.6. 

Figure 1 shows the convergence, shear, and both first and sec- 
ond flexion expected for a galatic-sized halo, assuming that the 
mass distribution is given by the Einasto, SIS, NEW, and Sersic 
profiles as a function of the angular separation from the lens cen- 
tre. An indicator of where the weak lensing effects are relevant 
is the Einstein radius of the SIS profile, which is 6e = 0.216 
for the halo that is being studied. As can be seen, the overall 
behaviours of the profiles are comparable, the differences be- 
tween the magnitudes of the lensing quantities being stronger at 
smaller angular separations, with the shear and second flexion 
showing a greater dissimilarity than the convergence and first 
flexion. These dissimilarities in the central region indicate that 
the Einasto profile may be differentiable from the other profiles 
in observational weak lensing studies. Our result agrees with 
Mamon et al. (2010), who compared the shear for the Einasto 
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Fig. 1. Comparison of the magnitude of convergence, shear, first and second flexion due to a dark matter halo of mass M200 = 
1 X lO'^/i"' Mq, assuming that the mass distribution is described by Einasto, NFW, and SIS profiles and the Sersic model. The 
lens and source planes are located at redshifts zl - 0.4 and zs = 0.92, respectively. The Einsten radius for the SIS is 6e - 0.216", 
indicating the angular distance when the weak lensing eff'ects are relevant. 



and NFW profiles finding that the major differences between the 
profiles are at small distances. As can be seen, the magnitudes of 
all properties at small distances are stronger for the SIS than the 
NFW profile and large distances the relation is inverted. In addi- 
tion, the magnitudes are stronger for the NFW profile than that 
of SIS, which is consistent with the slope for the SIS being -2, 
and for the NFW profile the inner slope being - 1 and the outer 
being -3. Interestingly, we also note that the lensing properties 
of the Einasto profile tend to be very similar to those of the NFW 
profile for increasing angular separation. 

In Figure 2, we plot dependence of the weak lensing properties 
on the concentration for the Einasto profile. We use the same 
mass as in the previous comparison, and values of c = 4, 8, 
12, 16, 20, 24 for the concentration. The concentration depen- 
dence of these properties is clearly non-linear, and very similar 
to that of the NFW profile investigated by Lasky & Fluke (2009). 
As with the NFW profile, the concentration dependence is far 
stronger at small separations. With increasing angular distance 
from the lensed image the dependence gets less pronounced 
and the curves become almost identical, except the convergence, 
which curves seem to be distinguishable at both small and large 
angular distances. The shear and second flexion seem to be more 
sensitive to the concentration variation than the convergence and 
first flexion, making the first two lensing properties useful tools 
for breaking the degeneracy between mass and concentration, if 



this occurs, as discussed by Lasky & Fluke (2009) for the NFW 
profile. 

We illustrate the index dependence of the Einasto profile in 
Figure 3. We again use the same halo mass M2(K) = 1 x 
10'^/;"' Mq, and vary the Einasto index, «, across the range 
1 < n < 9. It is evident, that the index dependence is greater 
at smaller angular separations for the lensing properties, becom- 
ing marginally distinguishable at increasing 6, with the only two 
exceptions being the exponential profile (n = 1), which can be 
clearly identified, and with a little bit more work the n = 2 pro- 
file can be distinguished as well. The convergence is the less de- 
pendent on « of all the lensing quantities; the shear and the first 
and second flexions depend more strongly on n than the conver- 
gence, near the lensed halo centre; this feature is very impor- 
tant, because it is near the lensed image that the lensing signal 
is stronger, making the first flexion, shear, and second flexion 
excellent tools for constraining n, whilst the convergence may 
be used to derive the halo mass. In contrast to the concentra- 
tion dependence of the Einasto profile, which is similar to that 
of the NFW profile, the index dependence of the Einasto profile 
differs from the corresponding dependence of the Sersic profile. 
Lasky & Fluke (2009) found that for the Sersic profile the index 
dependence is stronger in k and "F, and very weak in y and 0, 
whenever this dependence is present for the Sersic profile, its ef- 
fect decreases with boosted angular separations, where the lens- 
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ing signal is weaker, making it difficult to constraint the Sersic 
index. Thereby, we conclude that may be easier to constrain the 
index for an Einasto halo rather that for a Sersic halo. 

5. Summary and conclusions 

We have applied the Mellin transform technique to obtain 
closed-form expressions for the weak lensing properties of the 
Einasto profile. The expressions for the shear yix), first flex- 
ion 'F (x), and second flexion § (x) can be written in terms of 
the Fox H function for general values of the Einasto index n, 
and can simplified in terms of the Meijer G function for integer 
or half-integer values of n. We utilized the residue theorem to 
calculate specific power and logarithmic -power expansions for 
these expressions. The expansions permit us to study the asymp- 
totic behaviour of the weak lensing properties at small and large 
radii. Furthermore, we employed the Slater's theorem (Marichev 
1983) to derive an expression for the convergence /c (x) in terms 
of the generalized hypergeometric function, which is valid for 
half-integer values of n. This enables the other expressions for 
the lensing properties to be written in terms of the hypergeomet- 
ric function. 

We have examined in detail the convergence, shear, and first and 
second flexions for an Einasto profile and other profiles includ- 
ing the singular isothermal sphere, the NEW, and Sersic profiles. 
We found that the Einasto profile overall has a similar behaviour 
to these profiles. Nonetheless, this profile is clearly different 



from the others, particularly at small angular separations from 
the lens centre, where the lensing signal is stronger. At large 
angular separations, the Einasto profile behaves far very simi- 
larly to the NEW profile than the other profiles. We explored the 
dependence of the Einasto profile on the concentration parame- 
ter, our results indicating that it has a non-linear concentration 
dependence and that the shear and second flexion are more ef- 
fective indicators of the dependence than the convergence and 
the first flexion. In addition, we studied the Einasto index depen- 
dence. For this parameter, the dependence seems to be weaker 
than for the concentration, for which the shear, second, and first 
flexions seem to be more sensitive to the index dependence that 
the convergence. We note that the magnitude of the lensing prop- 
erties of a Sersic model are stronger than for the other profiles, 
this indicates that the profile selected to model the halo must be 
chosen with caution as discussed by Lasky & Fluke (2009). We 
note that the index dependence of an Einasto halo is stronger at 
small angular distances from the lens centre, which is the op- 
posite of the case for a Sersic halo for which at large angular 
distances the dependence is stronger This means that observa- 
tionally it is easier to constrain the value of the index for the 
Einasto profile rather than the Sersic model, because the lensing 
signal is stronger near the lensed image. 

The availability of analytical expressions for the Einasto-profile 
weak-lensing properties is of foremost importance, and consti- 
tutes an effort to foster the inclusion of this density profile in 
weak lensing modelling studies. There are several possible ap- 
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plications of our results for modelling studies. For example, one 
of them is the generation of the shear signal in weak lensing anal- 
yses, which can provide valuable information about the descrip- 
tion of the mass density profiles. This, in turn, places constraints 
on the model parameters such as mass, concentration, and par- 
ticularly the Einasto index, which is known to scale with mass 
and redshift according to A^-body simulations (Gao et al. 2008; 
Hayashi & White 2008), and for which our results could be used 
to verify this variation observationally. Likewise, weak flexion 
might be generated using our expressions and used to constrain 
the model parameters, particularly when halo substructure has to 
be proven, providing another scale within the haloes, where the 
behaviour of n can be studied. 

Here, we considered spherically symmetrical haloes, but haloes 
are far from ideal symmetric objects (see e.g. Shaw et al. (2006); 
Bettetal. (2007); Howell & Brainerd (2010)). Nevertheless, 
weak shear has proven successfully in studying deviations in the 
halo shape from spherical symmetry, such as the halo elliptic- 
ity (Hoekstra et al. 2004; Mandelbaum et al. 2006; Parker et al. 
2007; Evans & Bridle 2009; Oguri et al. 2010). Similarly, weak 
flexion has been proposed as a tool to investigate the halo el- 
lipticity (Hawken & Bridle 2009; Er & Schneider 201 1 ; Er et al. 
2011). Triaxiality is another aspect of the halo shape that has 
been explored using weak lensing (Oguri et al. 2005; Gavazzi 
2005; Sereno & Umetsu 2011; Feroz & Hobson 2012); ignor- 
ing the halo triaxiality can affect the parameter estimation in 
lens-rich clusters (Corless & King 2007; Corless et al. 2009), 



leading to the cluster appearing to be more massive and con- 
centrated, particularly, if its major axis is aligned with the 
line of sight (Hennawi et al. 2007; Oguri & Blandford 2009; 
Meneghetti et al. 2010). Some cluster studies where there are ap- 
parently lensing biases in the estimated concentration and mass 
include that of Broadhurst et al. (2008), who analysed the mass 
and concentration of four nearly relaxed clusters, that of the 
gravitational lens with the largest Einstein radius detected so 
fax- MACS J0717.5-I-3745 (Zitrin et al. 2009), and Oguri et al. 
(2009), who obtained the radial profile of four clusters com- 
bining lensing data from the Subaru telescope. The inclusion of 
our results in weak lensing ellipticity and triaxiality studies is 
straightforward, therefore it enables the possibility of investigat- 
ing the halo ellipticity and triaxiality with the Einasto profile. 
Our current knowledge of the structure of the Universe on large 
scales will be improved by new gravitational lensing surveys 
such as the Dark Energy Survey"* (DES), Euclid^, the Large 
Synoptic Survey Telescope (LSST)'', the James Webb Space 
Telescope^ KiDS^ Pan-STARRS'\ and WFIRST'". These sur- 



http ://w w w. darkenergy survey, org/ 
' http://sci.esa.int/euclid/ 
* http://www.lsst.org/lsst/ 
' http://www.jwst.nasa.gov/ 

http://kids.strw.leidenuniv.nl/ 
' http://pan-starrs.ifa.hawaii.edu/public/ 

http://wfirst.gsfc.nasa.gov/ 
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veys will provide more accurate measurements of weak lens- 
ing that could be modelled using the analytical expressions pre- 
sented in this work. 

This paper constitutes a further step in studying the properties of 
the Einasto profile using analytical means. In addition, it extends 
and complements the work of RM12, providing additional sim- 
plified expressions for their results. With this work, we hope to 
encourage the use of special functions such as the Fox H func- 
tion, the Meijer G function, and the generalized hypergeometric 
function in astronomy and astrophysics. 
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Appendix A: Simplified lialf-integer expressions of tlie Einasto profile lensing properties 

The eqs. (43)-(46) written in terms of Meijer G function can be reduced to expressions in terms of the generaUzed hypergeometric 
function (Luke 1969; Abramowitz & Stegun 1970) 



for half-integer values of n, for which the poles are simple, using the Slater's theorem (Marichev 1983) 
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where the asterisk in the gamma function indicates that the term k - j, which corresponds to F (0), must be replaced by 1 ; and in 
the hypergeometric function that the vector I + b^ - b must be reduced in size from qXo q- 1 . 
It is evident from eqs. (43) and (A.2) that the convergence may be written as 
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with bk the components of b given by eq. (43b). 

Inserting eq. (A. 3) into eqs. (32), (18) and (19), we obtain the expressions 
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Other important lensing quantities such as the cumulative surface mass density M {x), deflection angle a (x), and deflection potential 
ij/ (x) also can be written in terms of the generalized hypergeometric function 
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The eqs. (A.3)-(A.9) are ready for being use in numerical calculations because there are already several available numerical imple- 
mentations of the generalized hypergeometric function. Some software that includes an implementation of this function are the pro- 
prietary Maple, Mathematica, and Matlab, and freely available Sage, rapmath library, and the package hypergeo (Hankin & Lee 
2006) of the R language. 
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Appendix B: Series expansions of the shear and first and second flexions 

The expressions for the shear y {x) and first Tix) and second Q{x) flexions can be written as series expansions to study its asymptotic 
behaviour near zero. We apply the residue theorem to the contour integral in eq. (30) and obtain the explicit power or power- 
logarithmic series expansions depending on the multiplicity of the poles of the gamma functions \'{hj + BjS). Examples of specific 
applications of the residue theorem can be found in Baes & van Hese (2011), for the deprojected Sersic profile, and in RM12 for the 
projected Einasto profile. The general theorem for the Fox H function can be found in Kilbas & Saigo (1999). We encounter two 
cases: 

Case 1: if « is either non-rational or rational number pjq with an even denominator (and p, q are coprime), all poles are simple, so 
that the expansion takes the form of a power series 
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Case 2: if n is either an integer or rational number pjq with an odd denominator, some poles are of second order, then the expansion 
takes the form of a logarithmic -power series 
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with il/(k) the digamma function and ko 
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